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ABSTRACT:  Ocean  Acoustic  Waveguide  Remote  Sensing  (OAWRS)  has  recently  been  shown  to  be 
capable  of  instantaneously  imaging  and  continuously  monitoring  fish  populations  over  continental 
shelf-scale  areas,  covering  thousands  of  km2.  We  show  how  OAWRS  can  be  used  in  a  variety  of 
oceanic  ecosystems  to  remotely  assess  populations  and  study  the  behavior  of  fish  and  other  marine 
organisms,  such  as  Antarctic  krill,  to  help  the  study  of  marine  ecology  and  the  ecosystem-based 
approach  to  fisheries  management. 
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INTRODUCTION 

Fish  populations  often  comprise  the  largest  biomass 
component  in  a  productive  marine  ecosystem.  They 
typically  play  an  essential  role  in  inter-trophic  energy 
transport  (Cury  et  al.  2001),  and  serve  as  a  mainstay 
for  human  consumption  comprising  roughly  16%  of 
the  animal  protein  consumed  by  the  world's  popula¬ 
tion  (FAO  1997).  Despite  their  ecological  importance, 
there  is  substantial  evidence  that  fish  populations  are 
declining  worldwide  (Jackson  et  al.  2001,  Myers  & 
Worms  2003),  which  has  led  to  calls  for  an  eco¬ 
system  approach  to  fisheries  management  (Cury  et  al. 
2001,  Browman  &  Stergiou  2004,  Stergiou  &  Browman 
2005)  through  ecosystem-scale  sensing  and  monitor¬ 
ing  of  marine  habitats  (Godo  2009).  This  vision,  how¬ 
ever,  is  difficult  to  attain  with  conventional  methods 
(Simmonds  &  MacLennan  2005)  which  typically  rely 
on  in-situ  measurements  restricted  to  the  immediate 
vicinity  (tens  to  hundreds  of  meters)  of  slow  moving 
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research  vessels  and  greatly  undersample  the  ocean 
in  time  and  space  (Sund  1935,  Balls  1948,  Misund 
1997,  Medwin  &  Clay  1998,  Simmonds  &  MacLennan 
2005).  Recently,  a  method  known  as  Ocean  Acoustic 
Waveguide  Remote  Sensing  (OAWRS)  was  developed 
for  instantaneous  imaging  and  continuous  monitoring 
of  fish  populations  and  marine  life  over  continental 
shelf-scale  areas,  covering  thousands  of  km2,  at  an 
areal  rate  tens  of  thousands  to  millions  of  times 
greater  than  that  of  conventional  methods  (Makris  et 
al.  2006a).  Continuous  monitoring  with  OAWRS  pro¬ 
duces  unaliased  wide-area  movies  of  the  spatial  and 
temporal  distributions  of  fish  populations  that  can 
reveal  horizontal  behavioral  patterns  on  an  ecosystem 
scale.  This  may  enable  better  modeling  and  predic¬ 
tion  of  ecosystem  dynamics  as  well  as  correlation  with 
physical  and  biological  variables  including  those 
describing  oceanography,  climate,  food,  predation 
and  human  activity,  and  help  to  realize  ecosystem 
approaches  to  fisheries  management. 

©  Inter-Research  2009  •  www.int-res.com 


138 


Mar  Ecol  Prog  Ser  395:  137-160,  2009 


Ocean  waveguide  acoustics  has  been  used  for 
underwater  remote  sensing  in  all  oceans  for  on  the 
order  of  a  century  (Urick  1983).  It  takes  advantage  of 
the  fact  that  it  is  always  possible  to  form  an  acoustic 
waveguide  in  the  ocean  by  trapping  sound  between 
the  ocean- atmosphere  and  the  ocean-seabed  bound¬ 
aries,  as  is  described  in  numerous  textbooks  (e.g. 
Bergmann  1948,  Officer  1958,  Clay  &  Medwin  1977, 
Urick  1983,  Brekhovskikh  &  Lysanov  1991,  Frisk  1994, 
Jensen  et  al.  2000).  In  the  early  19th  century  the  prin¬ 
cipal  applications  were  naval  (Bergmann  1948,  Clay  & 
Medwin  1977,  Urick  1983).  The  world's  navies  still 
depend  upon  ocean  waveguide  acoustics  as  the  pri¬ 
mary  underwater  remote  sensing  tool  for  both  surface 
and  sub-surface  vessels  as  well  as  in  fixed  installations 
such  as  the  sound  surveillance  system  (SOSUS)  net¬ 
work  (Urick  1983,  Munk  et  al.  1995).  In  these  applica¬ 
tions,  horizontal  sensing  ranges  typically  span  dis¬ 
tances  many  orders  of  magnitude  larger  than  the  water 
depth  (Bergmann  1948,  Ewing  &  Worzel  1948,  Officer 
1958,  Tolstoy  &  Clay  1966,  Clay  &  Medwin  1977,  Urick 
1983,  Brekhovskikh  &  Lysanov  1991,  Frisk  1994, 
Jensen  et  al.  2000).  From  the  latter  part  of  the  19th  cen¬ 
tury  until  now,  ocean  waveguide  acoustics  has  been 
used  for  a  wide  variety  of  oceanographic  remote  sens¬ 
ing  applications,  including  quantification  of  ocean  and 
seabed  structure  and  passive  tracking  of  vocalizing 
marine  mammals  (Officer  1958,  Tolstoy  &  Clay  1966, 
Clay  &  Medwin  1977,  Urick  1983,  Brekhovskikh  & 
Lysanov  1991,  Frisk  1994,  Jensen  et  al.  2000,  Medwin 
2005). 

In  this  paper,  we  first  review  the  background  and 
some  of  the  recent  ecological  findings  made  with 
OAWRS.  We  describe  a  potential  method  for  remote 
species  classification,  which  we  demonstrate  with  field 
data.  We  then  show  how  OAWRS  may  be  used  in 
oceanic  ecosystems  to  remotely  assess  populations  and 
study  the  behavior  of  fish  and  other  pelagic  organisms, 
such  as  Antarctic  krill,  over  broad  temporal  and  spatial 
scales.  A  technical  description  of  the  OAWRS  approach 
appears  in  Appendices  A  to  F. 

BACKGROUND 

The  use  of  acoustics  to  detect  oceanic  fish  dates  back 
to  the  vertical  echo-sounders  introduced  in  the  1930s 
(Balls  1948,  1951).  During  the  1940s,  researchers  in  the 
University  of  California's  Division  of  War  Research 
(UCDWR)  noted  a  mid-water  layer  scattering  agent, 
which  was  later  called  'deep  scattering  layer'  (DSL) 
(Duvall  &  Christensen  1946,  Eyring  et  al.  1948,  Raitt 
1948).  After  bathypelagic  fish  with  gas-filled  swim- 
bladders  were  proposed  to  be  the  cause  of  this  DSL 
(Marshall  1951),  the  frequency  response  of  these  fish 


was  studied  to  identify  resonance  (Hersey  et  al.  1962, 
Andreeva  1964,  Andreeva  &  Chindonova  1964,  Mar¬ 
shall  &  Chapman  1964,  Andreeva  1965).  For  a  good 
historical  review  of  the  vast  literature  on  the  ability  of 
gas-filled  fish  bladders  to  scatter  sound  see  Hersey  & 
Backus  (1962)  and  Weston  (1967).  In  all  the  above 
experiments,  the  range  at  which  the  fish  were  imaged 
was  <500  m. 

Weston  &  Revie  (1971)  used  a  fixed  single-beam 
sonar  in  a  monostatic  setting  to  image  underwater 
sonar  returns  over  long  ranges  (>  10  km)  in  a  narrow 
horizontal  angular  sector.  They  observed  temporal 
variations  that  were  believed  to  be  consistent  with  fish 
migrations,  but  lacked  data  to  confirm  this.  Rusby  et  al. 
(1973)  generated  synthetic  aperture  images  of  the  con¬ 
tinental  shelf  environment  with  a  towed,  single-beam 
sidescan  sonar.  Each  synthetic  aperture  image 
required  hours  of  surveying,  which  led  to  high  spatio- 
temporal  undersampling  and  aliasing.  They  described 
features  as  possible  fish  groups  'only  when  the  shape 
of  the  groups  remain[ed]  sufficiently  distinctive  from 
run  to  run,'  which  would  bias  the  analysis  towards 
highly  static  population  distributions.  They  then 
guided  a  fishing  vessel  to  the  location  of  such  a  feature, 
where  the  vessel  made  a  large  fish  catch.  In  these  and 
other  earlier  long  range  experiments  (Nero  &  Huster 
1996,  Farmer  et  al.  1999),  independent  confirmation  of 
fish  was  not  available  by  simultaneous  measurements. 

Makris  et  al.  (2006,  2009)  used  a  horizontal  array  that 
formed  simultaneous  beams  over  a  360°  horizontal 
azimuth,  enabling  them  to  conduct  OAWRS  surveys  of 
marine  life  instantaneously  over  wide  areas,  tens  of 
thousands  of  km2.  With  regular  and  rapid  temporal 
image  updates,  they  were  able  to  work  in  a  true  Euler- 
ian  reference  frame  and  map  fish  distributions  without 
aliasing  in  space  or  time.  They  established  long- 
range  ocean  acoustics  as  a  method  for  detecting,  imag¬ 
ing,  and  estimating  fish  populations  (Appendix  B)  by 
coregistration  with  large  numbers  of  simultaneous 
OAWRS,  conventional  fish-finding  sonar  (CFFS)  and 
trawl  samples,  where  the  latter  provided  direct  species 
identification.  Large  numbers  of  simultaneous  mea¬ 
surements  are  necessary  for  confirmation  because  fish 
are  ubiquitous  in  continental  shelf  environments  and 
can  easily  be  found  accidentally  in  a  region  causing 
strong  acoustic  returns.  Non-simultaneous  correlations 
can  then  easily  be  spurious,  as  can  correlations  at  only 
a  single  or  very  small  number  of  spatial  locations.  This 
lesson  was  learned  with  geologic  features  of  the  sub¬ 
bottom,  which  are  also  ubiquitous  in  many  continental 
shelf  environments,  and  often  have  spurious  spatial 
correlation  with  acoustic  returns  caused  by  other 
mechanisms  (Ratilal  et  al.  2005). 

In  earlier  work  at  very  short  ranges  on  the  order  of 
the  water  depth,  -300  m,  and  so  with  conventional 
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direct-path  rather  than  waveguide  propagation  and 
sensing,  Isaacs  &  Schwartzlose  (1965)  used  a  U.S.  Navy 
mine  hunting  sonar  operating  on  the  southern  Califor¬ 
nia  continental  shelf  to  detect  strong  scatterers  which 
they  confirmed  to  be  fish  with  local  trawls. 

INVESTIGATING  FISH  SHOALING  BEHAVIOR 
OVER  ECOSYSTEM  SCALES 

The  OAWRS  approach  for  studying  marine  life  was 
first  demonstrated  in  2003  off  the  US  Continental  Shelf 
south  of  Long  Island,  NY,  i.e.  the  Mid-Atlantic  Bight. 
With  a  single  transmission  of  a  1  s  duration  linear  fre¬ 
quency  modulated  (LFM)  waveform,  OAWRS  sur¬ 
veyed  an  area  as  large  as  the  state  of  Connecticut  or 
New  Jersey  (Fig.  1).  The  imaging  is  effectively  instan¬ 
taneous  because  the  entire  region  is  surveyed  in  less 
time  than  it  takes  a  marine  organism  to  traverse  a 
single  OAWRS  resolution  cell  (Makris  et  al.  2006b). 
The  OAWRS  approach  was  used  again  with  the 
National  Marine  Fisheries  Service  (NMFS)  annual  her¬ 
ring  survey  of  the  Gulf  of  Maine  and  Georges  Bank  to 
study  herring  group  behavior  associated  with  spawn¬ 
ing  in  September-October  2006  (Makris  et  al.  2009a). 
In  both  the  OAWRS  experiments,  a  vertical  source 
array  transmitted  sound  in  the  frequency  range  of  390 
to  1400  Hz,  i.e.  near  swimbladder  resonance  for  many 
fish  species  in  the  survey  regions.  Echoes  scattered 
from  fish  were  received  by  a  towed  horizontal  receiv- 


Fig.  1.  Areal  coverage  of  a  single  Ocean  Acoustic  Wakeguide 
Remote  Sensing  (OAWRS)  transmission  during  the  2003  sur¬ 
vey  on  the  US  east  coast  continental  shelf.  An  area  of  60  km 
diameter  was  surveyed  every  40  s  (red  circle),  or  120  km 
every  80  s  (white  circle),  depending  on  ping  repetition  rate 
and  recording  time 
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Eastings  from  OAWRS  source  (km) 

Fig.  2.  Instantaneous  OAWRS  image  showing  fish  shoals  near 
the  continental  shelf  edge  100  km  south  of  Long  Island,  New 
York  (May  15,  2003,  10:36  Eastern  Daylight  Time  [EDT]). 
Dashed  white  lines  mark  depth  contours.  Receiver  array  reso¬ 
lution  decreases  as  viewing  directions  go  from  normal  (broad¬ 
side)  to  parallel  (endfire)  to  the  array  axis,  leading  to  blurring 
of  the  eastern  portion  of  the  northeastern  shoal.  Population 
density  estimation  employs  waveguide  propagation  and  scat¬ 
tering  models,  correction  for  OAWRS  areal  resolution,  and 
calibration  with  local  conventional  fish-finding  sonar  (CFFS) 
measurements  as  described  in  Appendices  A  to  F 


ing  array.  Instantaneous  snapshots  of  the  ocean  envi¬ 
ronment  over  thousands  of  km2  (Fig.  2)  were  then 
formed  by  charting  acoustic  returns  in  horizontal 
range  and  bearing  by  temporal  matched  filtering  and 
beamforming. 

With  the  'first  look'  of  OAWRS  in  the  Mid-Atlantic 
Bight  (Fig.  1)  it  was  possible  to  make  a  number  of  fun¬ 
damental  scientific  discoveries  about  the  (1)  instanta¬ 
neous  horizontal  structural  characteristics,  (2)  temporal 
evolution,  and  (3)  propagation  of  information  within 
very  large  fish  shoals  (Makris  et  al.  2006a).  These 
include  the  findings  that:  the  instantaneous  spatial  dis¬ 
tribution  of  fish  observed  follows  a  fractal  or  power  law 
process,  so  that  structural  similarity  exists  at  all  scales 
from  metres  to  tens  of  km  (previously  evidence  for 
structural  similarity  existed  only  for  small  scales 
<100  m;  Freon  &  Misund  1999);  large  shoals  (Fig.  2)  are 
far  more  horizontally  contiguous  in  2D  than  was  previ¬ 
ously  believed  based  on  ID  line  transect  methods 
which  sometimes  inaccurately  portray  them  as  disjoint 
clusters  (Pitcher  &  Parrish  1993,  Freon  &  Misund  1999); 
the  temporal  autocorrelation  scale  of  population 
change  within  a  very  large  shoal  is  extremely  short,  on 
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Fig.  3.  Areal  coverage  of  a  single  OAWRS  transmission  in  the  Gulf  of 
Maine,  2006).  A  region  of  100  km  diameter  (red  circle)  is  surveyed  every 
75  s.  Line  transects  of  National  Marine  Fisheries  Service  (NMFS) 
2  wk  survey  are  shown  in  yellow 


the  order  of  minutes,  which  is  why  fish  shoals  can  sud¬ 
denly  disappear  from  conventional  survey  vessels; 
temporal  fluctuations  in  shoal  population  also  follow  a 
power-law  process,  making  the  shoals  far  more  pre¬ 
dictable;  and  fish  density  waves  regularly  propagate 
information  over  km  scales,  3  orders  of  magnitude 
larger  than  previously  observed  (Shaw  1970,  Radakov 
1973),  at  speeds  an  order  of  magnitude  faster  than  fish 
can  swim,  which  apparently  help  large  shoals  remain 
cohesive.  These  observations  were  made  from  dis¬ 
tances  >10  km  from  the  shoals  with  sound  at  least  3 
orders  of  magnitude  less  intense  than  CFFS  due  to  the 
efficiency  of  ocean  acoustic  waveguide  propagation. 

General  predictions  about  animal  group  behavior 
believed  to  apply  in  nature  irrespective  of  species  (Vic- 
sek  et  al.  1995,  Toner  &  Tu  1998,  Reebs  2000,  Couzin  et 
al.  2005,  Buhl  et  al.  2006)  were  confirmed  by  monitor¬ 
ing  the  Georges  Bank  marine  ecosystem  (Fig.  3)  with 
OAWRS  in  the  fall  of  2006.  By  guantifying  the  forma¬ 
tion  process  of  vast  oceanic  herring  shoals  during 
spawning,  it  was  shown  that  (1)  a  rapid  transition  from 
disordered  to  highly  synchronized  behavior  occurs  as 
fish  population  density  reaches  a  critical  value;  (2) 
organized  group  migration  occurs  after  this  transition; 
and  (3)  small  sets  of  leaders  significantly  influence  the 
actions  of  much  larger  groups  (Makris  et  al.  2009a). 


The  spawning  process  was  found  to  follow  a 
regular  diurnal  pattern  in  space  and  time 
which  proved  to  be  difficult  to  detect  without 
continuous  wide-area  sensing  abilities  (Fig.  3). 
First,  pre-existing  populations  of  diffusely 
scattered  herring  reached  a  critical  density  at 
one  or  more  discrete  locations  near  the  north¬ 
ern  flank  of  Georges  Bank  just  before  sunset, 
apparently  in  response  to  diminishing  light 
level.  The  emergence  of  leading  clusters  then 
triggered  a  shoal-forming  convergence  wave 
(Fig.  4)  that  propagated  tens  of  km  in  tens  of 
min,  i.e.  at  speeds  an  order  of  magnitude 
greater  than  herring  can  swim  (Huse  &  Ona 
1996,  Fig.  3  of  Makris  et  al.  2009a).  Subse- 
guent  migrations  were  observed  towards 
southern  spawning  grounds  on  Georges  Bank, 
immediately  after  shoals  formed.  The  evi¬ 
dence  suggests  the  primary  biological  func¬ 
tion  of  the  shoals  is  a  prelude  to  synchronized 
spawning  and  the  shoals  form  in  deeper  water 
with  migrations  under  cover  of  darkness  to 
avoid  predators  (Makris  et  al.  2009a). 

As  Vicsek  et  al.  (1995)  theoretically  pre¬ 
dicted,  and  Buhl  et  al.  (2006)  showed  with  lab¬ 
oratory  experiments,  if  an  individual  assumes 
the  mean  speed  and  direction  of  those  in  its 
sphere  of  perception,  then  a  rapid  transition 
from  disordered  to  highly  synchronized 
behavior  occurs  when  a  critical  population  density  is 
attained.  This  may  be  understood  by  noting  that  as  the 
radii  of  perception  begin  to  overlap  with  increasing 
population  density,  chain  reactions  become  possible. 
Laboratory  experiments  (Reebs  2000)  and  simulations 
(Couzin  et  al.  2005)  showed  that  a  small  number  of 
individuals  can  significantly  influence  decision  mak¬ 
ing  in  large  groups. 

REMOTE  SENSING  OF  SWIMBLADDER  PROPERTIES 

For  fixed  fish  length  and  water  depth,  acoustic 
scattering  at,  or  near,  swimbladder  resonance  is  a 
strong  function  of  swimbladder  volume  (Fig.  5,  Supple¬ 
ment  1  at  www.int-res.com/articles/suppl/m395pl37_ 
app.pdf).  The  acoustic  scattering  can  vary  significantly 
across  species  and  so  may  be  helpful  in  species  classi¬ 
fication,  as  has  been  discussed  for  salmon  (Nero  & 
Huster  1996),  blue  whiting  (Love  1993),  Pacific  hake 
(Nero  et  al.  1998),  and  Atlantic  herring  (Nero  et  al. 
2004).  Swimbladder  volume  can  also  be  used  to  esti¬ 
mate  neutral  buoyancy  depth,  where  a  fish's  weight  is 
balanced  by  its  buoyancy  (Nero  et  al.  2004).  At  any 
given  depth,  neutral  buoyancy  reguires  the  swimblad¬ 
der  to  occupy  ~5%  of  the  total  fish  volume  (Jones  & 
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Fig.  4.  Example  of  regular  diur¬ 
nal  pattern  of  large  pre-spawn¬ 
ing  shoals  forming  near  sunset 
(18:08  h  EDT)  on  Oct  03,  2006. 
Top  panel:  1  hr  45  min  before 
sunset — no  shoal  present.  Bot¬ 
tom  panel:  10  min  after  sunset — 
large  shoal  present.  White 
dashed  lines  mark  depth  con¬ 
tours.  Population  density  esti¬ 
mation  employs  waveguide 
propagation  and  scattering 
models,  correction  for  OAWRS 
areal  resolution,  and  calibration 
with  local  CFFS  measurements 
as  described  in  Appendices  A  to 
F.  The  positive  vertical  axis 
points  16°  counter-clockwise  of 
true  north 


£ 

(U 

o 

c 

US 

to 

TD 

c= 

ea 

n 

1 

m 

w 

£ 


-30  -25  -20  -15  -10  -5  0  5  10  15 

Along-bank  distance  (km) 


Marshall  1953),  where  more  air  is  required  to  maintain 
this  ratio  as  depth  increases  due  to  the  compressive 
effects  of  increasing  pressure. 

We  used  OAWRS  to  determine  swimbladder  volume 
and  help  classify  the  species  of  shoaling  fish  observed 
during  the  2003  OAWRS  experiment  (Fig.  2)  (Makris 
et  al.  2006a).  The  mean  scattering  cross-section  of  a 
shoaling  fish  determined  with  OAWRS  2003  data  and 
local  CFFS  constraints  is  shown  in  Fig.  5  in  terms  of  tar¬ 
get  strength  (TS)  following  the  approach  given  in 
Appendix  E  for  3  frequencies  415,  925,  and  1325  Hz 
corresponding  to  the  center  frequencies  of  the  390  to 
440,  875  to  975  and  1250  to  1400  Hz  LFM  waveforms 
used  in  OAWRS  2003  (Makris  2003,  Makris  et  al. 
2006a).  The  least  squares  fit  between  the  empirically 
determined  OAWRS  TSs  and  those  determined  from  a 


standard  fish  scattering  model  (Love  model,  Supple¬ 
ment  1),  with  CFFS-measured  constraints  on  fish  habi¬ 
tation  depth  and  length,  leads  to  a  neutral  buoyancy 
depth  of  78  m  with  a  resonance  peak  at  -700  Hz  (black 
solid  line  in  Fig.  5),  consistent  with  the  CFFS-measured 
shoal  layer  at  70  to  90  m  depth.  While  buoyancy  depth 
need  not  correspond  to  habitation  depth,  the  corre¬ 
spondence  is  advantageous  because  it  enables  fish  to 
expend  minimal  energy  to  maintain  a  depth. 

While  the  ability  to  regulate  neutral  buoyancy  over  a 
wide  range  of  depths  is  usually  associated  with  physo- 
clist  fish  (closed  swimbladders),  physostome  fish  (swim- 
bladder  with  open  duct  to  the  intestine)  also  have  this 
ability  (Brawn  1962,  Sundnes  &  Sand  1975).  Physoclists 
are  known  to  slowly  regulate  the  swimbladder  volume 
through  the  blood  stream  (Jones  &  Marshall  1953)  and 


Fig.  5.  Target  strength  (TS)  data  for  3  frequencies  (cir¬ 
cles  ±  SD;  independent  samples;  415  Hz:  n  =  181,  SD  = 
0.7  dB;  925  Hz:  n  =  46,  SD  =  5.5  dB;  1325  Hz:  n  =  46,  SD 
=  5  dB;  see  Appendix  E)  corresponding  to  the  mean 
scattering  cross  section  of  a  shoaling  fish  species  in  the 
OAWRS  2003  experiment  constrained  by  local  CFFS. 
TS  frequency  curves  for  5  different  neutral  buoyancy 
depths  are  computed  with  Love's  model  (US  Navy  stan¬ 
dard  for  low  frequency  fish  scattering;  Supplement  1  at 
www.int-res . com/articles/suppl/m395p  1 3 7_app .pdf )  by 
depth-averaging  the  expected  scattering  cross  section 
of  an  individual  fish  over  the  layer  observed  by  CFFS 
for  the  Gaussian  length  distribution  (SD  15%  of  the 
28.6  cm  mean)  determined  by  CFFS.  The  least-squares 
best-fit  buoyancy  depth  between  measured  data  to 
Love-model  TS  is  given  by  the  black  solid  line.  If  the 
shoaling  fish  observed  in  the  70  to  90  m  layer  were  neu¬ 
trally  buoyant  closer  to  the  surface  (grey  line),  they 
would  scatter  far  too  weakly  below  1 .4  kHz  to  be  con¬ 
sistent  with  the  measured  OAWRS  and  CFFS  data 
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so  can  achieve  neutral  buoyancy  at  any  depth,  given 
enough  time  to  adjust  (Blaxter  &  Tytler  1978,  Harden 
Jones  &  Scholes  1985).  Physostomes,  on  the  other  hand, 
can  fill  their  swimbladders  by  gulping  air  at  the  surface 
(Blaxter  &  Batty  1984).  If  this  was  the  only  mechanism  for 
them  to  increase  the  volume  of  gas  in  their  swimbladder, 
then  their  maximum  neutral  buoyancy  depth  would  be 
limited  by  the  maximum  volume  they  could  gulp  at  the 
surface  (Thorne  &  Thomas  1990,  our  Fig.  SI).  Thorne  & 
Thomas  (1990)  measured  neutral  buoyancy  depths  of  up 
to  ~60  m  for  herring,  a  physostome,  and  suggest  addi¬ 
tional  gas  production  mechanisms,  such  as  from  bacteria 
in  the  digestive  tract  (Brawn  1962).  In  1  of  6  experiments, 
Brawn  (1962)  found  swimbladder  volume  build-up  by 
non-gulping  mechanisms  within  less  than  24  hr  in  her¬ 
ring  that  were  heavily  fed  with  copepods  containing  bac¬ 
teria.  Neutral  buoyancy  depths  of  40  to  50  m  have  been 
determined  from  near  resonance  scattering  data  (Nero  et 
al.  2004).  Char,  another  physostome,  have  also  been 
found  to  build  up  gas  in  their  swimbladder  even  when 
denied  access  to  the  surface  (Sundnes  &  Sand  1975).  The 
fish  imaged  by  OAWRS  were  consistently  observed  in  a 
70  to  90  m  layer  near  the  seafloor  for  >  12  h  by  CFFS.  The 
long  periods  that  the  fish  spent  at  these  depths  could  be 
sufficient  for  both  physoclists  and  physostomes  to  build 
up  gas  in  their  swimbladder. 

From  trawl  surveys  of  the  New  Jersey  continental  shelf 
taken  1  mo  prior  to  the  2003  experiment,  we  list  the  most 
probable  species  that  could  have  comprised  the  shoals 
imaged  by  OAWRS  based  on  freguency  of  catch:  At¬ 
lantic  herring,  scup,  hake,  black  sea  bass  and  ale  wife 
(NEFSC  Resource  Survey  Report  2003,  Makris  et  al. 
2006b).  While  dogfish  and  mackerel  are  also  found  in  the 
area  (NEFSC  Resource  Survey  Report  2003),  their  scat¬ 
tering  responses  are  expected  to  be  much  lower  than 
those  of  the  mentioned  fish  as  dogfish  and  mackerel  lack 
gas-filled  bladders.  Of  the  species  typically  found  in  the 
region,  only  Atlantic  herring  and  alewife  are  phy¬ 
sostomes,  while  scup,  hake  and  black  sea  bass  are 
physoclists.  Since  Atlantic  herring  is  the  most  common 
species  sampled  in  the  region  (NEFSC  Resource  Survey 
Report  2003),  and  the  only  one  known  to  form  such  large, 
extended  shoals,  they  are  most  likely  the  major  con¬ 
stituent  of  the  large  shoals  imaged  by  OAWRS  in  2003 
and  the  dominant  cause  of  scattering  measured  by  both 
OAWRS  and  CFFS.  The  fit  between  measured  TS  data 
and  the  Love  fish  scattering  model  constrains  the  sus¬ 
pected  herring  to  neutral  buoyancy  at  78  m  depth,  which 
suggests  that  they  should  have  a  mechanism  for  building 
up  gas  in  their  swimbladder,  or  they  have  damping 
mechanisms  which  lower  and  spread  the  resonance  and 
have  shallower  neutral  buoyancy  depth.  It  is  unlikely 
that  another,  less  abundant  species  could  have  domi¬ 
nated  the  scattering  measured  by  OAWRS  and  CFFS, 
since  this  would  require  'contaminants'  with  unrealisti¬ 


cally  large  target  strengths  to  follow  the  exact  spatial 
distributions  of  the  shoals  observed  at  both  OAWRS  and 
CFFS  frequencies.  Such  contamination  is  inconsistent 
with  trawl  data  (NEFSC  Resource  Survey  Report  2003), 
the  observed  preference  of  fish  to  shoal  among  similar¬ 
sized  individuals  (Pitcher  &  Parrish  1993),  and  stationar- 
ity  of  CFFS  and  OAWRS  scattering  measurements  (Ap¬ 
pendix  E),  which  show  fish  of  similar  length  and  target 
strength.  The  target  strengths  of  shoaling  fish  in  the 
OAWRS  2003  experiment  are  consistent  with  those 
measured  for  herring  in  the  2006  experiment  where 
trawl  samples  enabled  direct  species  identification 
(Makris  et  al.  2009a). 

POTENTIAL  ECOSYSTEM  EXPLORATION 

OAWRS  can  be  used  to  remotely  sense  a  variety  of 
fish  species  and  other  marine  organisms,  such  as 
Antarctic  krill,  over  broad  temporal  and  spatial  scales 
in  ocean  ecosystems  from  knowledge  of  (1)  the 
expected  scattering  cross  section  of  an  individual  (Sup¬ 
plement  1),  (2)  typical  population  densities,  (3)  acoustic 
propagation  (Appendix  D)  and  (4)  seafloor  scattering 
in  each  environment  (Appendix  F) .  Since  the  emphasis 
of  the  main  text  is  on  marine  biology  and  ecology,  we 
discuss  physical  acoustics  issues  in  Appendices  A  to  F 
and  Supplement  1. 

A  wide  range  of  transmission  frequencies  may  be 
used  for  OAWRS,  ranging  from  very  low  frequencies  of 
several  Hz  to  high  frequencies  of  tens  of  kHz.  Long- 
range  sound  propagation  in  the  ocean  is  less  attenuated 
at  lower  frequencies  and  is  less  sensitive  to  oceano¬ 
graphic  fluctuations.  Reverberation  from  the  seafloor 
tends  to  be  weaker  at  lower  frequencies. 

When  designing  an  OAWRS  system  we  need  to  con¬ 
sider  operating  at  frequencies  where  fish  scattering 
responses  are  high.  Optimal  OAWRS  frequencies 
should  be  low  enough  for  scattering  from  any  fish  to  be 
effectively  omni-directional,  to  make  OAWRS  insensi¬ 
tive  to  variations  in  fish  orientation,  which  is  typically 
not  the  case  in  CFFS.  Frequency  should  also  be  low 
enough  for  the  received  acoustic  field  scattered  from 
any  given  fish  to  be  expressible  as  the  product  of  the  3 
factors  'transmission  to  the  given  fish',  'scattering  from 
the  fish',  and  'transmission  from  the  fish'  (Appendix  C). 
At  CFFS  wavelengths,  this  factorization  is  typically  not 
possible  because  propagation  and  scattering  effects 
are  convolved  together  in  an  ocean  waveguide  (Ratilal 
et  al.  2002).  Frequencies  should  be  chosen  so  that 
acoustic  attenuation  from  propagation  through  the  fish 
is  negligible  even  over  long  ranges.  If  the  frequency  is 
too  low,  on  the  other  hand,  the  waveguide  may  no 
longer  support  modal  propagation  (Appendix  A),  and 
so  remote  sensing  may  become  inefficient. 
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Fig.  6.  Theragra  chalcogramma.  Modeled  target  strength  (TS) 
for  juvenile  and  adult  pollock  in  the  Gulf  of  Alaska  for  varying 
neutral  buoyancy  depths  (60  to  130  m).  The  expected  TS  at  a 
given  frequency  is  found  by  averaging  the  scattering  cross 
section  over  the  range  of  depths  and  mean  body  lengths  typi¬ 
cally  associated  with  these  fish  (Table  1).  The  body  lengths  of 
adults  and  juveniles  are  assumed  to  have  a  Gaussian  distribu¬ 
tion,  and  a  SD  of  10%  of  the  respective  means.  Details  of  the 
modeling  appear  in  Supplement  1 


Alaskan  pollock  Theragra  chalcogramma 

Alaskan  pollock  are  a  semi-pelagic  schooling  fish 
widely  distributed  in  the  North  Pacific  Ocean.  They 
comprise  the  largest  percentage  of  biomass  of  any  spe¬ 
cies  in  the  Bering  Sea1  and  are  an  important  predator 
that  feed  on  smaller  fish  and  also  cannibalize  their 
juveniles  (Dwyer  et  al.  1987).  So,  in  the  Gulf  of  Alaska, 
adults  are  more  often  found  near  the  bottom  (110  to 
130  m)  while  juveniles  prefer  to  stay  separate  above 
(70  to  100  m)  (Wilson  et  al.  2003).  We  use  the  parame¬ 
ters  summarized  in  Table  1  to  describe  the  pollock,  a 
physoclist  species  assumed  to  be  neutrally  buoyant  at 
the  depths  they  occupy. 

When  operating  at  the  resonance  peak  near  600  Hz 
(Fig.  6),  OAWRS  should  be  able  to  detect  adult  pollock 
schools  of  a  density  of  5  fish  nr2  with  a  signal  to  noise 
ratio  (SNR)  (Appendix  F)  of  30  to  40  dB  (Fig.  7). 
OAWRS  detections  of  pollock  should  then  span  a 
dynamic  range  of  30  to  40  dB  in  population  density, 
from  maximum  values  of  5  fish  nr2  (Table  1)  to  mini¬ 
mum  detectable  values  of  5  x  10-4  fish  m~2.  Single  pol¬ 
lock  should  then  be  observable  to  a  range  of  about 
3  km  with  the  OAWRS  2003  system.  Above  resonance, 
densities  of  at  least  10-2  fish  m“2  are  required,  but  this 
is  still  much  lower  than  the  typical  shoaling  densities  of 


^Alaska  Fisheries  Science  Center  NOAA  Fisheries.  Accessed 
26  Feb  2009.  www.afsc.noaa.gov/species/pollock.php 


Scattered  intensity  (dB) 
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pollock  (Table  1).  Even  at  these  higher  frequencies, 
detections  above  seafloor  scattering  should  span  a 
dynamic  range  of  at  least  2  to  3  orders  of  magnitude,  or 
20  to  30  dB,  in  population  density. 

Schools  of  juvenile  pollock  should  be  detectable  by 
OAWRS  at  resonance,  which  is  relatively  broad,  over  a 
dynamic  range  of  at  least  30  dB  (Fig.  7)  in  population 
density,  from  maximum  expected  values  of  10  fish  nr2 
(Table  1)  to  minimum  detectable  values  of  0.01  fish 


Peruvian  anchovy  Engraulis  ringens 

Peruvian  anchovy  play  an  essential  role  in  the  Hum- 
bolt  Current  upwelling  ecosystem  as  an  energy  inter¬ 
mediary  between  phytoplankton  and  large  predatory 
fish,  such  as  hake  and  horse  mackerel.  The  diets  and 
livelihood  of  seabirds,  marine  mammals,  livestock  and 
humans,  are  also  directly  tied  to  anchovy.  For  details 
on  the  Peruvian  anchovy  fishery  see  Niquen  &  Freon 
(2006). 

Anchovy  have  a  fairly  low  tolerance  to  temperature 
variations,  preferring  to  stay  in  the  13  to  23°C  range 
(Gutierrez  et  al.  2007).  Therefore,  their  geographic  and 
depth  distributions  are  highly  susceptible  to  climate 
change,  such  as  the  El  Nino  Southern  Oscillation 
(ENSO).  The  deepening  of  the  thermocline  and  the 
migration  of  food  sources  during  ENSO  periods  cause 
anchovies  to  descend  to  greater  depths  (Niquen  &  Diaz 
2002). The  average  depth  range  of  the  Peruvian  stock  is 


Range  (km) 

Fig.  7.  Theragra  chalcogramma.  Scattered  intensity  from 
adult  and  juvenile  pollock  and  the  seabottom  modeled  at 
600  Hz  (black  lines)  and  1  kHz  (grey  lines).  Pollock  are  as¬ 
sumed  to  be  uniformly  distributed  in  depth  layers  (adult:  110 
to  130  m;  juvenile:  70  to  100  m)  and  at  densities  of  5  adult  and 
10  juvenile  ind.  m-2  (Table  1).  A  description  of  the  modeling 
appears  in  Appendices  C,  D  and  F 


Fig.  8.  Engraulis  ringens.  Modeled  target  strength  for  an¬ 
chovy  off  the  Peruvian  coast  for  different  neutral  buoyancy 
depths.  Same  procedure  as  that  employed  in  Fig.  6.  Anchovy 
in  0  to  30  m  during  off-ENSO:  black  lines.  Anchovy  in  40  to 
70  m  during  ENSO:  grey  lines.  The  variation  in  neutral  buoy¬ 
ancy  corresponds  to  minor  axis  at  the  surface  expanding  by  a 
factor  of  ~1.4  to  2.2  (Fig.  SI  in  Supplement  1) 

0  to  30  m  (Bouchon  et  al.  2002,  Yanez  et  al.  2004) 
except  during  ENSO,  when  they  are  found  at  a  range 
of  40  to  70  m  (Blanco  et  al.  2002). 

Anchovy  have  an  average  length  of  14  to  16  cm 
(Bouchon  et  al.  2002)  (maximum  length  -20  cm; 
Gutierrez  et  al.  2007).  Catches  in  2005  showed  an  aver¬ 
age  length  of  14.5  cm  (Niquen  &  Freon  2006).  Anchovy 
are  known  to  form  huge  schools  (Tran  2003)  and  pack¬ 
ing  densities  of  115  (Graves  1977)  to  1312  (Smith  1970) 
fish  nr3  have  been  reported  for  Engraulis  mordox. 
Areal  densities  then  range  from  500  to  20  000  fish  nr2, 
given  typical  vertical  school  extents  of  4  (Holliday  & 
Larson  1979)  to  15  m  (Huntley  &  Zhou  2004).  E.  rin¬ 
gens,  E.  japonicus  and  E.  encrasocolus  display  similar 
schooling  characteristics  as  E.  mordox  (Huntley  &  Zhou 
2004). 

During  both  ENSO  and  off-ENSO  periods,  anchovy 
should  be  detectable  by  OAWRS  over  a  dynamic  range 
of  45  dB  in  population  density,  from  20  000  to  minimum 
detectable  values  of  0.6  fish  nr2,  which  is  much  lower 
than  typical  schooling  densities  (Table  1). 

When  operating  at  1.2  kHz,  where  uncertainties  in 
target  strengths  are  lowest  (Fig.  8),  OAWRS  should  be 
able  to  detect  anchovy  schools  with  an  areal  density  of 
20  000  fish  nr2  (15  m  thick  layer  of  fish  (Huntley  & 
Zhou  2004);  1  body-length  inter-fish  spacing  (Pitcher  & 
Partridge  1979);  Table  1)  with  an  SNR  of  45  dB  (Fig.  9) 
during  off-ENSO  periods.  During  ENSO,  OAWRS 
should  be  able  to  detect  anchovy  schools  of  20  000  fish 
nr2  with  an  SNR  of  45  dB  (Fig.  9)  when  operating  at 
2.4  kHz,  above  resonance  (Fig.  8). 
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Fig.  9.  Scattered  intensity  from  anchovy  during  off-ENSO 
(black  lines)  and  ENSO  periods  (grey  lines)  and  sea  bottom 
modeled  at  1.2  kHz  and  2.4  kHz  respectively,  and  for  a  maxi¬ 
mum  density  of  20  000  fish  m-2.  The  schools  are  assumed  to  be 
distributed  with  uniform  probability  within  the  top  30  m  (off- 
ENSO),  and  within  40  to  70  m  (ENSO)  depth 


Barents  Sea  capelin  Mallotus  villosus 

The  Barents  Sea  holds  the  largest  capelin  population 
in  the  world  and  this  species  is  the  largest  pelagic  fish 
component  of  the  region,  with  a  biomass  reaching  6  to 
8  x  106  t  (Gjosaeter  1998).  Capelin  play  a  key  role  in 
energy  conversion  from  zooplankton  to  higher  trophic 
predators,  such  as  haddock,  harp  seal,  Northeast  Arc¬ 
tic  cod,  whales,  and  sea  birds  (Gjosaeter  1998),  and  is 
also  important  to  human  and  livestock  consumption 
(Froese  &  Pauly  2009).  Capelin  undergo  drastic  diel 
vertical  migration  patterns,  occupying  depths  of  30  to 
60  m  at  night  and  125  to  175  m  at  day  (Jorgensen 
2004).  Packing  densities  of  capelin  in  the  Barents  Sea 
vary  from  (day)  1.4  to  15  and  (night)  0.8  to  3.5  fish  nr3 
(Serberov  1985).  Using  vertical  school  extents  of  (day) 
50  and  (night)  30  m  (Jorgensen  2004),  we  obtain  mean 
school  densities  of  (day)  70  to  750  and  (night)  24  to  105 
fish  nT2.  Capelin  body  lengths  average  16  cm  (Jor¬ 
gensen  2003,  Gauthier  &  Horne  2004)  (maximum:  20 
cm  for  females,  25  cm  for  males;  Froese  &  Pauly  2009). 
Capelin  are  characterized  by  a  shorter  swimbladder 
than  most  fish  with  a  major  axis  that  is  approximately 
15%  of  their  length  (Jorgensen  2003,  Gauthier  & 
Horne  2004). 

When  operating  at  the  resonance  peak  near  1.2  kHz 
(Fig.  10),  OAWRS  should  be  able  to  detect  capelin 
schools  of  a  density  of  105  fish  nr2  (night;  Table  1) 
with  an  SNR  of  30  dB  (Fig.  11).  Capelin  should  then 
be  detectable  by  OAWRS  over  a  dynamic  range  of 
30  dB  in  population  density  between  105  and  0.1  fish 


Frequency  (Hz) 

Fig.  10.  Mallotus  villosus.  Modeled  target  strength  (TS)  for 
capelin  at  2  different  water  layers  (Table  1)  in  the  Barents 
Sea,  for  different  neutral  buoyancy  depths.  Shallow  (night): 
black  lines;  deep  (day):  grey  lines.  Same  procedure  as  that 
employed  in  Fig.  6 

m“2,  which  is  much  lower  than  typical  schooling  den¬ 
sities.  Daytime  schools  of  capelin  (750  fish  nr2,  Table  1) 
should  be  detectable  by  OAWRS  at  resonance 
(2.4  kHz)  with  an  SNR  of  at  least  45  dB  (Fig.  11).  Day¬ 
time  schools  of  capelin  should  then  be  detectable  by 
OAWRS  over  a  dynamic  range  of  45  dB  in  population 
density  from  750  to  a  minimum  value  of  0.02  fish  m-2, 
which  are  much  lower  than  typical  schooling  densities 
(Table  1). 


Range  (km) 

Fig.  11.  Mallotus  villosus.  Scattered  intensity  from  capelin  and 
seabottom  modeled  at  1.2  kHz  (night:  black  lines)  and  at 
2.4  kHz  (day:  grey  lines).  Capelin  are  assumed  to  have  pack¬ 
ing  densities  of  (night)  105  and  (day)  750  fish  m-2,  and  to  be 
uniformly  distributed  in  a  depth  layer  (night:  30  to  60  m;  day: 
125  to  174  ) 
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Southern  blue  whiting  Micromesistius  australis 

The  continental  shelf  and  adjacent  waters  off  South¬ 
ern  Patagonia  constitute  one  of  the  main  fishing 
grounds  of  the  Argentine  Sea.  The  Southern  blue  whit¬ 
ing  is  the  most  abundant  demersal-mesopelagic  fish  in 
this  region  (Cassia  2000),  and  is  the  third  largest  spe¬ 
cies  caught  by  local  fishing  fleets  (Cassia  2000,  Macchi 
&  Pajaro  2005),  after  Argentine  hake  and  hoki.  South¬ 
ern  blue  whiting  typically  occupy  water  depths  of  100 
to  200  m  (Agnew  2002)  along  the  Patagonia  continen¬ 
tal  slope  and  shelf  break  and  provide  a  major  food 
source  for  larger  higher-trophic  predators,  such  as 
Argentine  hake,  whales,  sea  lions,  as  well  as  sea-birds, 
including  the  yellow-eyed  penguin  and  black-browed 
albatross  (Thompson  et  al.  1998,  Cherel  et  al.  1999, 
Jackson  et  al.  2000). 

The  TS  of  southern  blue  whiting  is  modeled  by 
assuming  a  uniform  vertical  distribution  over  a  50  m 
thick  layer  from  150  to  200  m  depth  (Cassia  2000). 
Since  whitings  are  physoclists,  they  are  expected  to  be 
neutrally  buoyant  over  a  similar  depth  range.  Trawl 
samples  from  typical  spawning  stocks  (Cassia  2000) 
show  a  bimodal  distribution  of  lengths  associated  with 
age  and  sex  (mean  length;  male:  39  cm;  female:  51  cm 
(Cassia  2000).  The  areal  density  of  1  fish  nr2  (Table  1) 
is  the  average  density  derived  from  trawl  catches  over 
a  large  area  covering  hundreds  of  km2  (O'Driscoll  et  al. 
2005).  Southern  blue  whiting  are  known  to  mass  in 
dense  aggregations  wherein  the  school  densities  may 
be  much  higher  than  those  reported  in  Table  1. 

When  operating  at  the  resonance  peak  near  800  Hz 
(Fig.  12,  black  lines),  OAWRS  should  be  able  to  detect 
female-dominated  schools  of  an  areal  density  of  1  fish 


Fig.  12.  Micromesistius  australis.  Modeled  target  strength 
(TS)  for  southern  blue  whiting  in  the  Argentine-Falklands  re¬ 
gion  for  different  neutral  buoyancy  depths.  Black  lines:  female; 
grey  lines:  males.  Same  procedure  as  that  employed  in  Fig.  6 


m“2  with  an  SNR  of  30  dB  (Fig.  13).  Female-dominated 
schools  should  then  be  detectable  by  OAWRS  over  a 
dynamic  range  of  at  least  30  dB  in  population  density, 
from  values  of  1  to  0.001  fish  nr2,  which  is  much  lower 
than  typical  schooling  densities  (Table  1).  Male- 
dominated  shoals  of  typical  density  1  fish  nr2  should 
be  detectable  by  OAWRS  at  resonance  (~1  kHz;  Fig  12, 
grey  lines)  with  an  SNR  of  25  dB  (Fig.  13).  Male-domi¬ 
nated  schools  should  then  be  detectable  by  OAWRS 
over  a  dynamic  range  of  at  least  25  dB  in  population 
density  from  1  to  a  minimum  detectable  value  of 
0.003  fish  nr2  which  is  much  lower  than  typical  school¬ 
ing  densities  (Table  1).  It  should  then  be  possible  to 
observe  both  single  males  and  females  up  to  a  range  of 
about  2  km. 


Argentine  hake  Merluccius  hubbsi 

Argentine  hake  is  a  mid-trophic  level  species  provid¬ 
ing  forage  for  whales,  sea  lions,  penguins,  sea-birds, 
and  other  fauna  native  to  the  Patagonian  ecosystem 
(Alonso  et  al.  2000).  Argentine  hake  are  known  to  feed 
on  smaller  fish  such  as  anchovies,  juvenile  hake, 
southern  blue  whiting,  squids,  and  macrozooplankton 
(FAO  2009).  Argentine  hake  are  typically  found  in  the 
coastal  and  continental  shelf  environments  of 
Falkland- Malvinas  Islands  between  the  100  to  200  m 
isobaths.  Fishing  pressures  on  Argentine  hake  in  the 
1990s  caused  shifts  to  deeper  water  and  lower  density 
spawning  populations  (Macchi  et  al.  2007). 

Argentine  hake  are  known  to  amass  in  large  shoals 
in  coastal  waters  at  depths  between  60  and  90  m. 


5  10  15  20  25  30 


Range  (km) 

Fig.  13.  Micromesistius  australis.  Scattered  intensity  from 
whiting  and  seabottom  modeled  at  800  Hz  (female:  black 
lines)  and  at  1  kHz  (male:  grey  lines).  The  whiting  are  as¬ 
sumed  to  have  a  packing  density  of  1  fish  m-2,  and  to  be  dis¬ 
tributed  from  150  to  200  m  depth  with  uniform  probability 
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Spawning  occurs  from  October  to  February,  with  peak 
activity  in  January.  Early  in  the  spawning  season,  hake 
populations  are  dominated  by  juveniles  of  both  sexes 
with  a  mean  length  ~28  cm  (Macchi  et  al.  2007). 
Towards  the  end  of  the  spawning  season,  larger  adult 
hake  (mean  lengths  ~38  cm)  accumulate  along  with 
the  juveniles  (Macchi  et  al.  2007).  The  densities 
reported  in  Table  1  are  average  densities  derived  from 
trawl  catches  over  a  large  area  covering  4  x  104  km2 
(Macchi  et  al.  2007).  Hake  are  known  to  mass  in  dense 
aggregations  wherein  the  school  densities  may  be 
much  higher  than  those  reported  in  Table  1.  Both  juve¬ 
nile  and  adult  hake  are  physoclist,  and  are  assumed  to 
be  neutrally  buoyant  at  60  to  90  m  depth. 

When  operating  at  the  resonance  peak  near  650  Hz 
(Fig.  14,  black  lines),  OAWRS  should  be  able  to  image 
adult  hake  shoals  of  population  density  of  0.5  fish  nr  2 
(Table  1)  with  an  SNR  of  25  dB  (Fig.  15).  OAWRS 
detections  of  adult  hake  should  then  span  a  dynamic 
range  of  at  least  25  dB  in  population  density  from  0.5 
(Table  1)  to  a  minimum  detectable  value  of  0.002  fish 
m-2.  Above  resonance,  densities  of  at  least  0.02  fish  m-2 
are  reguired,  but  this  is  still  much  lower  than  the 
average  shoaling  densities  (Table  1),  so  that  OAWRS 
imagery  should  span  a  dynamic  range  of  at  least  2 
orders  of  magnitude,  or  20  dB,  in  population  density. 

Schools  of  juvenile  hake  with  a  population  density  of 
0.3  fish  nr2  should  be  detectable  by  OAWRS  at  their 
resonance  of  900  Hz  (grey  lines  in  Fig.  14)  with  an  SNR 
of  20  dB  (Fig.  15).  OAWRS  detections  of  juvenile  hake 
should  then  span  a  dynamic  range  of  at  least  20  dB  in 
population  density,  from  0.3  fish  nr  2  to  a  minimum 
detectable  value  of  0.003  fish  nr2. 


Fig.  14.  Merluccius  hubbsi.  Modeled  target  strength  (TS)  for 
hake  in  the  beginning  (corresponding  to  a  majority  of  juve¬ 
niles:  grey  lines)  and  end  (corresponding  to  a  majority  of 
adults:  black  lines)  of  the  spawning  season,  off  Argentina. 

Same  procedure  as  that  employed  in  Fig.  6 


Fig.  15.  Merluccius  hubbsi.  Scattered  intensity  from  hake  and 
seabottom  modeled  at  650  Hz  (adults:  black  lines)  and  900  Hz 
(juveniles:  grey  lines).  During  the  day,  hake  are  assumed  to 
be  uniformly  distributed  in  depth  layers  (adults:  60  to  90  m, 
packing  density  0.5  fish  m-2;  juveniles:  60  to  90  m,  density 
0.3  fish  m-2) 

Atlantic  bluefin  tuna  Thunnus  thynnus 

Atlantic  bluefin  tuna  are  large  top -predators  that 
feed  on  mid-sized  fish.  Their  size  and  speed  allows 
them  to  evade  most  predators,  with  the  exception  of 
sharks,  larger  toothed  whales  and  humans  (Cascorbi  & 
Chabot  2004).  As  physoclists,  neutral  buoyancy  is 
expected  at  their  common  swimming  depths  of  10  to 
40  m  (Lutcavage  2000).  A  typical  tuna  school  is  ~20  m 
in  diameter  (Newlands  et  al.  2006),  with  an  inter-fish 
separation  of  1  body  length  (areal  density  of  0.25  fish 
m-2)  (Partridge  et  al.  1983). 

For  the  OAWRS  2003  system,  where  operating  fre- 
guencies  are  above  the  resonance  at  ~50  Hz  (Fig.  16), 


Fig.  16.  Thunnus  thynnus.  Modeled  target  strength  (TS)  for 
bluefin  tuna  with  different  neutral  buoyancy  depths  (10  to 
40  m).  Same  procedure  as  that  employed  in  Fig.  6 
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Fig.  17.  Thunnus  thynnus.  Scattered  intensity  from  bluefin 
tuna  and  seabottom  modeled  at  950  Hz.  Tuna  packing  density 
within  an  OAWRS  resolution  footprint  is  assumed  to  be  range- 
dependent  and  is  given  by  the  ratio  'area  occupied  by  the 
tuna  school'  :  'area  of  the  resolution  footprint'  x  'school  den¬ 
sity  of  0.25  fish  m-2'.  For  example,  for  the  OAWRS  2003  sys¬ 
tem,  the  effective  density  of  a  school  within  the  OAWRS  foot¬ 
print  would  be  0.025  fish  m-2  at  10  km.  At  20  km,  the  footprint 
area  doubles  and  the  effective  density  reduces  to  0.0125  fish 
m-2.  This  is  reflected  in  the  figure  as  a  reduction  in  the  signal 
to  noise  ratio  as  the  range  increases.  The  tuna  are  assumed  to 
be  uniformly  distributed  in  a  layer  from  0  to  30  m  water  depth 

typical  tuna  schools  occupy  areas  smaller  than  an 
OAWRS  resolution  cell  for  ranges  >70  m.  The  effective 
areal  density  of  a  single  school  within  an  OAWRS 
resolution  footprint  is  then  given  by  the  ratio:  'area 
occupied  by  the  tuna  school'  :  'area  of  the  footprint)'  x 
'school  density  of  0.25  fish  m-2'.  For  example,  at  a 
range  of  10  km  the  effective  areal  density  is  -0.025  fish 
m“2.  Since  the  effective  density  is  a  function  of  the 
OAWRS  resolution  footprint  area,  the  SNR  varies  with 
range  (Fig.  17)  from  15  to  20  dB. 

An  OAWRS  system  could  be  designed  to  operate 
near  resonance  (50  Hz)  and  with  sufficiently  high  reso¬ 
lution  so  that  the  'area  occupied  by  a  single  typical 
tuna  school'  >  'OAWRS  resolution  cell'  (Appendix  C). 
In  this  scenario,  OAWRS  could  image  typical  tuna 
schools  with  an  SNR  of  50  dB,  given  the  higher  TS  and 
effective  density.  OAWRS  detections  of  tuna  should 
then  span  a  dynamic  range  of  50  dB  in  population  den¬ 
sity,  from  0.25  fish  m“2  (Table  1)  to  minimum  detectable 
values  of  3  x  10-6  fish  m“2,  or  -1  fish  in  a  600  x  600  m 
area. 

The  wide-area  spatial  coverage  and  continuous  tem¬ 
poral  monitoring  of  OAWRS  can  be  an  asset  in  study¬ 
ing  the  behavioral  dynamics  and  spatial  distributions 
of  fast- swimming,  highly  migratory  pelagic  fish.  Since 
bluefin  tuna  swim  at  speeds  of  4  km  hr1  (Itoh  et  al. 
2003)  and  occupy  broad  geographic  scales,  they  are 
difficult  to  survey  with  conventional  methods  such  as 


electronic  tagging,  satellite,  or  spotter  plane  aerial 
imaging  which  are  limited  to  studying  individual  fish 
or  surficial  populations  (Lutcavage  2000,  Newlands  et 
al.  2006). 

The  ability  of  OAWRS  to  image  and  continuously 
track  small,  but  rapidly  swimming  schools  of  marine 
creatures  of  -100  m  extent  was  demonstrated  in  2003 
(Makris  et  al.  2006a, b).  For  example,  a  small  compact 
group,  located  at  9  km  south  and  3  km  east  of  the 
OAWRS  source,  was  observed  traveling  north  at  -5  km 
hr1,  consistent  with  a  typical  tuna  school  (Fig.  18).  Sim¬ 
ilar  aggregations  were  observed  within  a  1.5  km  radius 
exhibiting  the  morphological,  dynamical  spatial  dis¬ 
tributions,  speeds  and  uniform  trajectories  typical  of 
fast-swimming  tuna  schools.  The  OAWRS-measured 
densities  for  these  groups  of  -5  x  10-4  fish  nr2,  after 
compensating  for  the  expected  TS  of  an  individual  at 
415  Hz  ( — 17  dB),  follows  expectations  for  tuna. 

Antarctic  krill  Euphausia  superba 

Antarctic  krill  is  widely  distributed  in  the  Antarctic 
ocean  (Nicol  &  Endo  1997)  and  plays  a  key  role  as  the 
primary  source  of  sustenance  for  many  species  of 
whales,  seals,  birds,  fish  and  sguid  (Miller  &  Hampton 
1989). 

It  may  be  possible  for  OAWRS  to  detect  and  image 
large  swarms  of  krill,  given  a  properly  designed 
OAWRS  system  (Chen  2008).  Krill  differs  only  slightly 
in  compressibility  and  density  from  the  surrounding 
seawater  as  it  does  not  have  any  gas-filled  cavities. 
This  makes  krill  extremely  weak  sound  scatterers  at 
OAWRS  frequencies  (Fig.  19)  compared  to  swimblad- 
dered  fish.  Therefore,  krill  concentrations  need  to  be 
much  higher  than  those  of  swimbladdered  fish  for  suc¬ 
cessful  OAWRS  imaging.  Fortunately,  typical  swarm 
densities  of  Antarctic  krill  exceed  103  ind.  nr3,  and  can 
reach  up  to  105  ind.  m“3  in  superswarms  (Hamner  & 
Hamner  2000).  While  these  high  densities  help  to  off¬ 
set  the  weak  scattering  of  an  individual,  effective  use 
of  OAWRS  to  image  krill  will  likely  require  taking 
advantage  of  the  vertical  sound  speed  structure  of  the 
polar  environments.  This  sound  speed  structure  makes 
it  possible  to  design  an  OAWRS  system  that  channels 
sound  in  well  known  and  widely  used  refract-surface- 
reflected  (RSR)  paths  that  do  not  interact  with  the 
seafloor,  but  only  with  the  ocean-atmosphere  inter¬ 
face  to  form  a  waveguide  (Urick  1983;  Supplement  1). 
Since  krill  are  mainly  found  in  the  upper  water  column 
(<100  m)  (Hamner  et  al.  1983,  Brierley  &  Watkins  2000, 
Brierley  et  al.  2002),  RSR  paths  can  insonify  the  krill 
well  with  minimal  bottom  interference.  Long  range 
sensing  systems  have  used  purely  waterborne  refrac¬ 
tive  paths  for  remote  sensing  of  submerged  objects  and 
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Fig.  19.  Euphausia  superba.  Target  strength  (TS)  for  different 
sizes  of  krill.  The  TS  is  calculated  by  Rayleigh-Born  approxi¬ 
mation  (Morse  &  Ingard  1968)  given  the  geometry  and 
acoustic  properties  of  krill  (Lawson  et  al.  2006) 


seamounts  in  scenarios  more  challenging  than  those 
proposed  here.  For  example,  during  a  field  experiment 
at  the  Mid- Atlantic  Ridge,  narrow  sound  beams  were 
guided  by  refraction  over  ranges  of  >100  km,  to  almost 
graze  the  ocean  bottom  and  sea  surface  at  ray  vertex, 
pass  through  narrow  ridges  and  accurately  impinge  on 
the  preordained  geological  features  under  investiga¬ 
tion  (Makris  et  al.  1995,  ONR  1996,  Smith  &  Cushman 
1997,  Harding  et  al.  1998,  Preston  2000,  Medwin  2005). 

The  scattered  field  from  krill  is  compared  with 
seafloor  scattering  at  10  kHz  in  Fig.  20  for  both  deep 
ocean  (2000  m)  and  continental  shelf  environments 
(200  m). 

An  OAWRS  system  operating  at  10  kHz  should  be 
able  to  detect  swarms  of  2  cm  krill  at  a  range  of  5  km 
when  their  packing  density  is  >0.1  ind.  m“3  for  deep 
and  >1  ind.  nr3  for  shallow  waters.  Since  typical 
swarming  densities  are  much  higher,  OAWRS  should 
be  able  to  image  typical  krill  swarms  with  a  dynamic 
range  of  at  least  40  dB  for  deep  ocean,  and  30  dB  for 
continental  shelf  environments  (Fig.  20). 

For  smaller  krill  of  1  cm  length,  OAWRS  imaging  is 
expected  to  be  less  favorable  in  continental- shelf  envi¬ 
ronments,  unless  densities  are  at  least  105  ind.  nr3 
(superswarms,  see  above).  A  transmitter  array  with 
lower  sidelobe  levels  may  resolve  this  problem  even 
for  typical  densities  of  1000  ind.  nr3. 


Summary 

The  dynamic  range  expected  in  OAWRS  imaging  is 
summarized  in  Fig.  21  for  a  variety  of  ecologically  sig¬ 
nificant  fish  species,  and  in  Fig.  22  for  Antarctic  krill. 
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Deep  water  environment;  Depth  =  2000  m 


Shallow  water  environment;  Depth  =  200  m 


Fig.  20.  Euphausia  superba.  Scattered  intensity  from  krill  and 
seabottom  for  (a)  deep  ocean  and  (b)  continental  shelf  envi¬ 
ronments.  Krill  individuals  are  assumed  to  have  lengths  from 
1  to  4  cm,  a  typical  packing  density  of  1000  ind.  m-3  (Hamner 
&  Hamner  2000),  and  a  uniform  distribution  from  0  to  100  m 
depth.  Black  and  grey  lines  correspond  to  scattered  intensity 
computed  using  the  winter  and  summer  profiles  in  Supple¬ 
ment  1  Fig.  S2,  respectively 

For  all  fish  species  examined,  typical  shoaling  densi¬ 
ties  (Fig.  21)  are  at  least  2  orders  of  magnitude  greater 
than  the  minimum  densities  detectable  by  OAWRS, 
making  them  viable  candidates  for  future  wide-area 
surveys.  For  Antarctic  krill  of  at  least  2  cm  length, 
OAWRS  should  be  capable  of  imaging  typical  swarm¬ 
ing  densities  (Fig.  22)  with  a  dynamic  range  at  least  2 
orders  of  magnitude.  With  appropriate  system  design, 
it  may  be  possible  to  image  smaller  krill  with  OAWRS. 


CONCLUSIONS 

We  reviewed  a  number  of  recent  findings  in  marine 
ecology  related  to  the  behavior  of  vast  oceanic  fish 
shoals  that  were  made  possible  by  Ocean  Acoustic 
Waveguide  Remote  Sensing  (OAWRS),  a  technigue 


Fig.  21.  Comparison  of  dynamic  ranges  in  population  density 
expected  in  wide-area  surveys  for  the  8  fish  species  dis¬ 
cussed.  Lower  end  of  vertical  bar  =  expected  minimum  indi¬ 
vidual  densities  observable  with  OAWRS.  Upper  end  =  maxi¬ 
mum  areal  fish  densities  from  historical  observations.  Shaded: 
typical  shoaling  density  from  historical  observations 


Fig.  22.  Expected  dynamic  range  in  population  density  for 
Antarctic  krill  as  a  function  of  length  in  wide-area  surveys. 
Lower  end  of  vertical  bar  =  expected  minimum  individual 
density  observable  with  OAWRS  from  modeling  and  simula¬ 
tion  (Chen  2008).  Upper  end  =  maximum  densities  from  his¬ 
torical  observations.  Shaded  =  typical  swarming  densities 
from  historical  observations 
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capable  of  instantaneously  imaging  and  continuously 
monitoring  fish  populations  over  continental  shelf- scale 
areas,  spanning  thousands  of  km2.  OAWRS  can  be 
used  in  oceanic  ecosystems  to  remotely  assess  popula¬ 
tions  and  study  the  behavior  of  fish  and  other  marine 
organisms  such  as  Antarctic  krill.  The  approaches  pre¬ 
sented  here  for  wide-area  and  continuous-time  moni¬ 
toring  of  pelagic  species  may  help  to  meet  the  signifi¬ 
cant  demands  of  the  ecosystem-based  approach  to 
research  and  conservation  in  marine  biology  (Garcia  & 
Cochrane  2001,  Pikitch  et  al.  2004,  Browman  &  Ster- 
giou  2004,  Stergiou  &  Browman  2005). 

Currently,  OAWRS  has  been  deployed  from  moving 
research  vessels.  In  the  future,  it  will  likely  also  be 
deployed  at  fixed  locations  to  enable  continuous  long 
term  monitoring  of  oceanic  ecosystems  and  their  varia¬ 
tions.  A  precedent  already  exists  for  this  in  the  atmos¬ 
phere  where  fixed  Doppler-weather  radar  (DWR)  sta¬ 
tions  have  been  used  to  instantaneously  image  bird 
populations  and  study  population  distributions,  migra¬ 
tions  and  behavior  over  wide  areas.  Continuous 
DWR  monitoring  has  been  instrumental  in  avian  con¬ 
servation  and  ecosystem-based  resource  management 
(Gauthreaux  &  Belser  2003,  2005,  Ruth  et  al.  2005). 
With  freguent  use,  OAWRS  can  likewise  play  a  signifi¬ 
cant  role  in  scientific  exploration,  ecosystem  manage¬ 
ment  and  conservation  in  the  world's  oceans. 
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Appendix  A.  Ocean  acoustic  waveguides 


Remote  sensing  with  ocean  waveguide  acoustics  takes  ad¬ 
vantage  of  the  fact  that  it  is  always  possible  to  form  an 
acoustic  waveguide  in  the  ocean  by  trapping  sound  between 
the  ocean-atmosphere  and  the  ocean-seabed  boundaries 
(Bergmann  1948,  Ewing  &  Worzel  1948,  Officer  1958,  Tol¬ 
stoy  &  Clay  1966,  Clay  &  Medwin  1977,  Urick  1983, 
Brekhovskikh  &  Lysanov  1991,  Frisk  1994,  Jensen  et  al.  2000). 

In  continental  shelf  environments,  most  applications  of 
long  range  ocean  acoustics  involve  waveguides  that  are 
formed  by  interaction  of  sound  with  both  of  these  bound¬ 
aries  (Officer  1958,  Chapter  3;  Jensen  et  al.  2000,  Chapter  1; 
Clay  &  Medwin  1977,  Chapter  9;  Urick  1983,  Chapter  6; 
Brekhovskikh  &  Lysanov  1991,  Chapter  5).  Even  though  this 
is  often  referred  to  as  'shallow-water  acoustics'  or  'littoral 
acoustics'  the  acoustic  wavelength  is  much  smaller  than  the 
ocean  depth  so  that  acoustic  transmission  involves  many  prop¬ 
agating  modes  (Fig.  Al).  This  type  of  boundary-interacting, 
continental-shelf  propagation  is  considered  for  various  fish 
types  in  'Potential  ecosystem  exploration'  in  the  main  paper. 

Sometimes,  it  is  possible  to  form  acoustic  waveguides  in 
the  ocean  where  sound  does  not  interact  with  the  ocean- 
atmosphere  boundary,  ocean-seabed  boundary  or  both,  by 
appropriate  experimental  design  involving  restriction  to 
sufficiently  shallow  propagation  angles  (Officer  1958, 
Tolstoy  &  Clay  1966,  Clay  &  Medwin  1977,  Urick  1983, 
Brekhovskikh  &  Lysanov  1991,  Frisk  1994,  Jensen  et  al. 
2000).  In  the  deep  ocean  in  mid-latitudes,  for  example,  the 
combination  of  increasing  pressure  with  depth  and  rela¬ 
tively  warm  water  near  the  surface  typically  leads  to  a 
sound  speed  minimum  at  -1000  m  depth  where  the  acoustic 
field  can  be  trapped  and  propagate  thousands  of  km  with¬ 


out  interacting  with  either  boundary.  In  polar  environments, 
cold  water  near  the  surface  leads  to  sound  speed  minima  at, 
or  near,  the  surface;  so  a  waveguide  can  be  formed  by 
refraction  at  depth  and  surface  reflection  with  sound  rays 
that  never  interact  with  the  seafloor  boundary.  Such  sound 
speed  structures  have  been  exploited  over  the  years  in 
many  remote  sensing  applications  (Bergman  1948,  Officer 
1958,  Clay  &  Medwin  1977,  Urick  1983,  Brekhovskikh  & 
Lysanov  1991,  Frisk  1994,  Jensen  et  al.  2000).  In  our  analy¬ 
sis  for  krill  in  'Potential  ecosystem  exploration',  we  employ 
the  latter  approach,  known  as  refract-surface-reflect  (RSR; 
Urick  1983)  propagation  in  an  experiment  designed  for 
polar  environments  (see  Supplement  1  at  www.int-res.com/ 
articles/suppl/m395p  137_app.pdf ) . 

Regardless  of  the  specific  ocean  waveguide,  there  are 
many  standard  approaches  for  determining  the  acoustic 
field  in  an  ocean  waveguide,  including  those  using  normal 
mode,  wave  number  integration,  parabolic  equation  and  ray 
tracing  (Jensen  et  al.  2000)  formulations.  Similarly,  there  are 
many  standard  models  using  each  of  these  approaches,  e.g. 
Kraken,  Ocean  Accoustics  and  Seismics  Exploration  Syn¬ 
thesis  (OASES),  Range-dependent  Acoustic  Model  (RAM) 
and  Generic  Sonar  Model  (GSM)  (Porter  &  Reiss  1985, 
Jensen  et  al.  2000,  Collins  1993),  which  have  been  carefully 
benchmarked  over  the  years  (Buckingham  &  Tolstoy  1990, 
Collins  1990,  Jensen  &  Ferla  1990,  Thomson  1990,  Thomson 
et  al.  1990,  Stephen  1990,  Westwood  1990).  In  'Potential 
ecosystem  exploration',  we  use  the  US  Navy  standard  RAM 
parabolic  equation  model  to  determine  transmission  loss  in 
all  the  continental  shelf  environments  investigated  in  the 
present  study,  as  described  in  Appendix  D. 


Sea-surface  Pressure  =  0 


Sea-surface 


Fig.  Al.  Sketch  illustrating  bound¬ 
ary-interacting,  long-range,  modal 
propagation  in  a  typical  conti¬ 
nental  shelf  waveguide.  Spheri¬ 
cally  spreading  waves  from  a  point 
source  are  multiply  reflected  from 
the  ocean-atmosphere  and  ocean- 
seafloor  boundaries  to  form  vertical 
modes  that  propagate  horizontally. 
In  an  iso-sound  speed  layer,  each 
mode  can  be  expressed  as  a  vertical 
standing  wave  formed  by  the  inter¬ 
ference  of  an  up  and  down-going 
plane  wave  of  fixed  horizontal  graz¬ 
ing  angle  determined  by  the  layer's 
boundary  conditions.  The  sketch 
shows  modes  and  equivalent  plane 
waves  for  a  canonical  iso-sound 
speed  continental  shelf  environ¬ 
ment,  known  as  a  Pekeris  wave¬ 
guide 
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Appendix  B.  Experimental  and  theoretical  statistics  of  OAWRS  intensity  images 


In  the  Ocean  Acoustic  Waveguide  Remote  Sensing 
(OAWRS)  experiments  of  2003  and  2006  (Makris  et  al. 
2006a, b,  2009a, b)  acoustic  returns  measured  by  OAWRS 
were  found  to  obey  circular  complex  Gaussian  random 
(CCGR;  Goodman  1985)  field  fluctuations  following  theory 
(Bergmann  1948,  Dyer  1970,  Makris  1996,  Makris  et  al. 
2006a, b,  2009a, b)  and  many  previous  ocean  acoustic  exper¬ 
iments  (Dyson  et  al.  1976,  Lynch  et  al.  2003,  Fredricks  et  al. 
2005,  Ratilal  et  al.  2005,  Andrews  et  al.  2009b,  Tang  et  al. 
2008).  The  instantaneous  intensity  I  of  a  CCGR  field  fol¬ 
lows  the  exponential  distribution,  while  averaged  intensity 
(Goodman  1985)  and  the  log  of  averaged  intensity  (Makris 
1996)  follow  the  gamma  and  exponential-gamma  distribu¬ 
tions,  respectively,  with  first  and  second  moments  that  can 
be  analytically  expressed  in  terms  of  sample  size  p  and 
expected  intensity  (I)  (Makris  1996).  The  standard  deviation 
(SD)  of  the  log  of  averaged  intensity  from  a  CCGR  field  is 

a  =  (101og10e)VC(2,ia)  (Bl) 

where  £  is  Riemann's  zeta  function,  and  is  approximately 
given  by  o  =  4.34/l/p,  for  p  >  3  (Makris  1996)  so  that  sta¬ 
tionary  averaging  over  p  independent  samples  asymptoti¬ 
cally  reduces  the  SD  by  >/iTr  a  fact  that  is  exploited  in  most 
imaging  systems  which  employ  CCGR  fields,  including 
OAWRS,  to  obtain  low  variance  images. 

In  the  OAWRS  2003  experiment,  the  empirically  deter¬ 
mined  SD  of  the  log  of  averaged  intensity  from  OAWRS 


imagery  is  on  the  order  of  1  dB  at  a  given  OAWRS  pixel.  For 
example,  a  ~  1.3  dB  SD  is  measured  for  the  log  of  averaged 
intensity  at  any  pixel,  as  shown  in  Fig.  Bla  where  no  large 
fish  groups  are  present  and  the  scattered  returns  leading  to 
the  observed  intensity  measurements  are  consistent  with 
returns  from  the  seafloor.  The  same  ~1.3  dB  SD  is  also  mea¬ 
sured  for  OAWRS  imaging  of  a  large  fish  shoal  (Fig.  Blb,c). 
In  both  cases,  the  low  SD  are  the  result  of  10  sample  (5  tem¬ 
poral,  2  range)  averaging,  and  closely  follow  the  1.4  dB  the¬ 
oretically  expected  from  Eq.  (Bl),  as  has  been  noted  by 
Makris  et  al.  (2006a, b).  The  -20  dB  contrast  between  fish 
shoals  and  background  seafloor  scattering  (Fig.  Bl)  is  much 
larger  than  the  -1.3  dB  SD,  which  is  low  enough  to  make 
details  in  shoal  morphology  discernable.  Log-transformed 
intensity  is  used  to  present  OAWRS  imagery  data  since  pat¬ 
tern  recognition  in  intensity  images  formed  from  CCGR 
field  data  is  optimized  by  matching  in  the  log  domain  where 
noise  is  signal-independent  and  the  variance  is  stabilized 
(Makris  1995)  by  homomorphic  transformation.  Similar  SD 
for  the  log  of  measured  intensity  in  the  OAWRS  2006  exper¬ 
iment  were  obtained  by  stationary  averaging  (Makris  et  al. 
2009b). 

In  the  OAWRS  2003  and  2006  experiments,  averaged 
intensity  returns  from  the  seafloor  exhibited  a  trend  of 
smooth  decay  with  range  (Fig.  Bl,  Fig.  S5  in  Makris  et  al. 
2009b)  following  theoretical  expectations  for  uniformly  dis¬ 
tributed  waveguide  scatter ers  (Bucker  1970,  Ellis  1995, 
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Fig.  Bl.  (a)  Averaged  intensity  measured  by  OAWRS  in  the 
absence  of  prominent  fish  shoals,  normalized  to  unit  source 
power,  from  the  OAWRS  2003  survey  shows  a  trend  of 
smooth  decay  with  range.  Returns  are  consistent  with  scat¬ 
tering  from  the  seafloor.  Error  bars  =  experimentally  deter¬ 
mined  SD  of  -1.3  dB  for  the  standard  10  sample  (5-ping  and 
2-range-cell)  intensity  average  employed  in  OAWRS  2003. 
(b)  Instantaneous  OAWRS  image  of  fish  population  density 
showing  a  large  fish  shoal  (10:25  h  Eastern  daylight  time 
[EDT];  May  14,  2003).  Solid  black  line  =  transect  through 
the  shoal  along  which  averaged  intensity  is  shown  in  (a). 
This  image  is  typical  of  thousands  collected  during  the  2003 
OAWRS  survey  (Makris  2003).  (c)  Averaged  intensity  mea¬ 
sured  by  OAWRS  along  the  transect  through  a  large  fish 
shoal  in  (b),  with  experimentally  determined  SD  of  -1.3  dB 
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Appendix  B  (continued) 


Makris  &  Ratilal  2001)  and  many  measurements  of  seafloor 
reverberation  and  transmission  in  the  ocean  (Bucker  1970, 
Makris  et  al.  1995,  Smith  &  Cushman  1997,  Harding  et  al. 
1998,  Hines  et  al.  1998,  Preston  2000,  Ratilal  et  al.  2005, 
Andrews  et  al.  2009b,  Makris  et  al.  2009b).  The  smooth 
trend  is  expected  because  the  many  independent  modal 
contributions  from  random  propagation  and  scattering  that 
lead  to  CCGR  field  statistics  by  the  central  limit  theorem 
(CLT)  also  lead  to  a  lack  of  coherent  modal  interference 
structure  (Bucker  1970,  Ellis  1995,  Makris  &  Ratilal  2001, 
Chen  et  al.  2005,  Ratilal  &  Makris  2005,  Makris  et  al.  2006b, 
Chen  et  al.  2008).  Randomization  of  the  acoustic  field  arises 
from  both  fluctuations  in  the  water  column  due  to  diverse 
phenomena  such  as  internal  waves,  eddies,  turbulence, 
boundary  roughness  at  the  sea  surface  and  seafloor,  as  well 


as  fluctuations  in  the  properties  of  scatterers  such  as  posi¬ 
tion,  aspect  and  composition. 

The  instantaneous  intensity  of  conventional  fish-finding 
sonar  (CFFS)  returns  from  fish  vary  significantly  as  a  func¬ 
tion  of  fish  aspect  and  also  follow  CCGR  statistics  (Dahl  & 
Mathisen  1983)  as  a  consequence  of  the  CLT  (Clay  &  Med- 
win  1977).  It  has  also  been  noted  (Makris  1996)  that  as  a 
consequence  of  the  CLT,  the  same  CCGR  field,  averaged 
intensity  and  log  averaged  intensity  statistics  are  found  in 
CFFS  scattering  from  fish  and  in  ocean-acoustic  waveguide 
transmission  scintillation,  where  stationary  averaging  is  typ¬ 
ically  needed  in  all  cases  to  produce  low  variance  images. 
The  5.6  dB  instantaneous  intensity  SD,  for  example,  is  often 
too  high  for  many  imaging  applications  and  is  typically 
reduced  by  stationary  averaging  (Makris  1996). 


Appendix  C.  Expected  intensity  in  OAWRS  imaging  of  fish  groups 


The  expected  square  magnitude  of  the  field,  (lcj)s(pc)l2)r  proportional  to  instantaneous  intensity,  scattered  from  AT indepen¬ 
dent  and  identically  distributed  fish  with  random  position,  orientation  and  scattering  properties  within  the  OAWRS  resolution 
footprint  of  area  A(pc),  centered  at  horizontal  location  pc  can  be  expressed  as  (Andrews  et  al.  2009a) 


101og10(l<|>s(pc)l2}  =  SL  +  TLA(pc,  z0,H)  +  101og10 


+  101og10  (na) 


(Cl) 


where  SL  is  the  source  level  normalization,  TLA  is  a  transmission  loss  area  term  describing  the  expected  second  moment  of 
depth  averaged  propagation  to  and  from  the  fish  layer  integrated  over  the  resolution  footprint  of  the  OAWRS  system,  S(f)  is 
the  random  scatter  function  of  a  fish  in  the  group,  k  is  the  acoustic  wavenumber,  and  (na)  =  N/A( pc)  is  the  expected  areal  fish 
density  within  the  spatially  varying  resolution  footprint. 

The  third  term  on  the  right  hand  side  of  Eq.  (Cl)  is  defined  as  the  target  strength  (TS)  corresponding  to  the  expected  scat¬ 
tering  cross  section  of  a  fish  in  the  group.  The  TLA  term,  a  function  of  center  depth  z0  and  thickness  of  the  fish  layer  H,  can 
be  expressed  as 


TLA(pc,  z0,H)  =  101og10  %(pc,z0,H) 


X(Pc,Zo  ,H)  = 

zs=z0+H/2 

J^(pc)  J  JJ  (4?t)4  I  G(r  Ips.Zsi/.cfr^j.dfr^jGfps.Zglroi^cjr^J.dlr^))  l2P(c(rw),d(rlv))P(ps,zs)dzsd2psdc(rlv)dd(rw) 

A(pc)  zs~zo~H/2  c(rw),d(rw) 

where  G(r\ps,zs]f,c(rw),d{rw))  and  G(ps,zs  I r0-,f,c(rw),d(rw))  are  Green  functions  describing  random  waveguide  propagation  to 
and  from  the  fish,  P(c(rw),d{rw))  is  the  joint  probability  distribution  of  sound  speed  c  and  seawater  density  d  in  the  water  col¬ 
umn  at  any  point  rwin  the  propagation  path  and  P(ps,zs)  is  the  probability  of  finding  a  fish  at  ps,zs.  For  a  uniform  distribution 

1  1 

of  fish  within  the  OAWRS  resolution  footprint  P(ps,zs)  = - ,  so  that 

H  A(pc) 

^  zs=zo+H/2 

X(Pc'zo'H)  =  J  —  J  (47c)4(lG(rlpSlzs;/1c(rw),d(rJ)G(pSlzslr0;/,c(rM,)1d(rJ)l2|ps,zs)dzsd2ps  (C3) 

^(Pc)  zs-zo-^/2 

where  the  conditional  expectation  is  over  the  water  column  sound  speed  and  density  random  variables. 


Appendix  D.  Transmission  loss  over  the  OAWRS  resolution  footprint 


The  transmission  loss  area  (TLA)  term  in  Eq.  (Cl)  is 
computed  using  the  US  Navy  standard  range-dependent 
acoustic  model  (RAM)  (Collins  1993)  in  conjunction  with 
measured  oceanographic  data  including  sound  speed  of  the 
water  column  and  sediment,  sediment  density  and  attenua¬ 
tion,  and  bathymetry. 

The  TLA  computed  for  all  the  environments  considered  in 
this  paper  are  displayed  in  Fig.  Dla,b  and  show  the  trend  of 
smooth  decay  with  range  expected  from  theory  (Bucker 
1970,  Ellis  1995,  Makris  &  Ratilal  2001,  Chen  et  al.  2005), 
measured  seafloor  reverberation  (Fig.  Bla,  Galinde  et  al. 


2008)  and  past  experiments  (Lynch  et  al.  2003,  Fredricks  et 
al.  2005,  Ratilal  et  al.  2005,  Tang  et  al.  2008).  They  also 
exhibit  low  SDs  of  ~1  to  1.5  dB  over  expected  fish  shoal 
depths  (Table  1). 

To  determine  TLA,  the  Green  functions  in  Eq.  (C3) 
are  computed  with  the  RAM  parabolic  equation  model 
(Collins  1993,  Jensen  et  al.  2000)  for  an  ocean-acoustic 
waveguide.  The  conditional  expectation  value  in  Eq.  (C3) 
is  determined  by  averaging  100  Monte-Carlo  realizations, 
where  the  Green  functions  are  computed  along  the  entire 
propagation  path  in  range  and  depth  for  each  realization. 
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Appendix  D  (continued) 


Each  Monte-Carlo  realization  itself  employs  a  different 
sound  speed  depth  profile  every  500  m  (Chen  et  al.  2005, 
Andrews  et  al.  2009b)  along  the  propagation  path.  Mea¬ 
sured  sound  speed  profiles  from  the  OAWRS  2003  experi¬ 
ment  (Fig.  Die)  are  used  to  generate  TLA  in  Fig.  Dla,  and 
from  the  Argo  Database  (2009)  (Fig.  Did)  to  generate  TLA 
in  Fig.  Dl(b).  The  SD  of  TLA  over  a  depth  layer  H,  is 
defined  as 


J  z=z0+H/2 

a  (TLA)  =  J—  J  (TLA0(z)-TLA(H,z0))2dz  (Dl) 

V  ^  z=z0-H/2 

where 

TLA0(z)  =  101og10  J  (4n)“  (D2) 

A(Pc) 

(lG(rlps,z;/,c(rw)4(rJ)G(lpSlzlroi/,c(rJ,d(rJ)l2|ps,z)d2ps 


Fig.  Dl.  (a)  Expected  2-way  transmission  loss  area  term  (TLA)  computed  by  Monte-Carlo  simulation  for  the  OAWRS  2003 
experimental  environment  (continental  shelf  south  of  Long  Island,  NY),  for  the  H  =  30  m  and  z0  =  65  m  where  fish  shoals  were 
observed  by  CFFS.  Error  bars  show  the  TLA  SD  (o)  of  1  to  1.5  dB  over  the  depths  of  this  fish  shoal  layer,  (b)  Expected  2-way 
TLA  for  the  environments  considered  in  'Potential  ecosystem  exploration'  in  the  main  paper.  For  each  environment,  the  mean 
TLA  is  computed  by  averaging  over  the  expected  fish  shoal  depth  (Table  1).  Error  bars  indicate  the  SD  in  TLA  over  these 
depths.  The  water  depths  in  the  different  environments  are  (1)  150  m  for  Bering  Sea,  Gulf  of  Mexico,  Peru,  Barents  Sea,  (2) 
180  m  for  Gulf  of  Maine,  and  (3)  200  m  for  Argentina,  Antarctica,  (c)  Sound  speed  profiles  measured  during  the  OAWRS  2003 
survey  in  the  continental  shelf  south  of  Long  Island,  NY,  and  used  to  compute  TLA  in  (a),  (d)  Measured  sound  speed  profiles 
for  various  continental  shelf  environments  (Argo  Database  2009)  used  to  compute  TLA  in  (b).  The  environments  span  the 
canonical  cases  of  upward  refracting  (Antarctic),  well-mixed  (Gulf  of  Mexico),  downward  refracting  (Peru),  and  a  deep  water 

sound  speed  minimum  (Gulf  of  Maine) 
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Appendix  E.  Empirical  estimation  of  target  strength  and  areal  population  density  from  OAWRS  data 


The  target  strength  (TS)  at  OAWRS  frequencies,  TSqawrsf 
is  estimated  using  (1)  OAWRS  measurements  of  averaged 
intensity,  (2)  modeled  TLA,  and  (3)  simultaneous  CFFS  mea¬ 
surements  of  number  density.  Stationary  averaging  is  then 
employed  over  hundreds  of  OAWRS  and  CFFS  samples  to 
obtain  a  low  variance  TSqawrs  estimate,  using  OAWRS  fish 
population  density  maps  such  as  those  shown  in  Figs.  2  & 
Bib. 

The  maximum  likelihood  estimate  (MLE) 


TSoAWRs(psr-f)  =  101og10 


101og10  l<Mps,/)l2  - 


(El) 


101og10  %(pc,z0,H,  f)  -  101og10  na(pc)  -  SL 


is  obtained  from  invariance  of  the  MLE  (Kay  1993)  as  pre¬ 
scribed  by  evaluation  of  parameters  at  their  corresponding 
MLEs  in  Eq.  (Cl).  While  the  first  2  terms  on  the  right  hand 
side  of  Eq.  (El)  are  estimated  from  OAWRS  measurements 
of  averaged  intensity  (Appendix  B)  and  modeled  TLA 
(Appendix  D),  respectively,  the  areal  fish  number  density  na 
is  estimated  from  simultaneous  CFFS  measurements  as 


101og10  na(pc)  =  SSc(pc)-TSc  M 

where  SSC  is  CFFS-measured  scattering  strength  and  TSC  is 
the  target  strength  of  an  individual  fish  at  the  CFFS  operat¬ 
ing  frequency. 


Table  El.  SD  of  different  terms  on  the  right  hand  side  of 
Eqs.  El  and  E2  before  and  after  stationary  averaging  over 
OAWRS  pixels  (n  =  JL81).  Note  that  SD  of  the  CFFS  target 
strength  estimate  TSC  is  not  affected  by  averaging  over 
OAWRS  data 


Terms  in  Eqs.  (El)  &  (E2) 

SD  per 

OAWRS  pixel 
(dB) 

SD  after 
stationary 
averaging 
(dB) 

101og10UM'^7)l2 

1.3 

0.09 

W}og10x{pc,z0,H,f) 

0.5  to  1 

0.07 

SSc  (Pc) 

1.5 

0.11 

TSC 

0.65 

0.65 

TSoawrs 

2.3 

0.67 

We  now  explain  how  the  SD  of  the  fish  target  strength 
estimate  at  415  Hz  from  OAWRS  2003  data  shown  in  Fig.  5 
in  the  main  paper  was  determined  and  how  it  is  consistent 
with  theory  and  other  measurements.  At  a  given  OAWRS 
2003  pixel,  the  variance  (i.e.  SD2)  of  TSoawrs  is  the  sum 
of  the  variances  of  each  term  on  the  right  hand  side  of 
Eq.  (El),  given  their  independence.  For  the  first  2  terms, 
101og10  lcj)s(pc,/)  l2and  101og10£(p c,z0,H,f),  the  SD  are  1.3  dB 
(Appendix  B)  and  1  dB  (Appendix  D),  respectively.  The  vari¬ 
ance  of  the  third  term,  101og10na(pc) ,  is  the  sum  of  the  vari¬ 
ances  of  SSC  and  TSC ,  given  the  independence  of  the  terms 
in  Eq.  (E2).  The  SD  of  SSC  and  TSC  are  1.5  dB  and  0.65  dB 
per  OAWRS  pixel,  respectively,  from  CFFS  measurements 
made  during  the  OAWRS  2003  experiment.  The  resulting 
SD  of  TSoawrs,  per  OAWRS  pixel,  is  then  2.3  dB. 

Estimates  of  TSqawrs  based  on  many  independent  OAWRS 
pixels  within  regions  of  statistically  stationary  fish  popula¬ 
tions  are  then  averaged  to  reduce  the  SD.  During  the 
OAWRS  2003  experiment,  both  OAWRS  and  CFFS  co-regis- 
tered  many  fish  shoals  with  statistically  stationary  popula¬ 
tions,  such  as  the  one  shown  in  Fig.  Bl.  The  SD  of  all  the 
terms  on  the  right  hand  side  of  Eqs.  (El)  &  (E2)  after  station¬ 
ary  averaging  over  181  OAWRS  samples,  obtained  within 
such  stationary  populations,  are  summarized  in  Table  El. 
Note  that  the  SD  of  TSC  cannot  be  reduced  by  averaging 
OAWRS  data  because  it  is  assumed  not  to  vary  across 
OAWRS  samples.  After  stationary  averagingjhe  theoretical 
and  empirical  SD  of  0.67  dB  (Table  El)  for  TSoawrs  agrees 
well  with  the  empirically  determined  value  of  0.7  dB  shown 
in  Table  E2  &  Fig.  5. 

To  estimate  OAWRS  areal  number  densities  over  the  wide 
areas  shown  in  Figs.  3  in  the  main  paper  &  Bib,  Eq.  (Cl) 
is  again  employed  by  now  grouping  together  all  terms 
except  101og10na(pc),  andjusing  the  empirically  estimated 
TSqawrs-  Extrapolating  TSoawrs  in  an  OAWRS  image  to 
spatial  locations  where  there  are  no  CFFS  measurements  is 
valid  when  the  fish  obey  stationary  random  processes  in 
their  spatial  distribution  and  scattering  properties. 

The  variance  of  the  OAWRS  areal  number  density,  at 
any  giveryjDAWRS  pixel,  is  the  sum  of  the  variances  of 
101og10  l(|)s(pc,/)  I2  10 log  10 %(pc ,z0,H,f)  and  TSoawrs.  Using 
the  SD  in  Tables  El  &  E2,  OAWRS  2003  estimates  of  number 
density  then  have  a  SD  of  1 .5  dB  for  statistically  stationary 
fish  populations. 


Table  E2.  Empirically  estimated  target  strength  (TS)  at  415  Hz  for  OAWRS  2003  experiment.  Six  transects  through  statistically 
stationary  fish  populations  co-registered  by  OAWRS  and  CFFS  are  used  to  compute  the  least  squares  estimate  (effectively 
same  as  MLE  for  the  given  data)  and  standard  deviation  of  TSqawrs-  nv  =  number  density.  EDT  =  Eastern  Daylight  Time 


Day 

(May  2003) 

Time 

(EDT) 

Transect 
length  through 
shoal  (m) 

Average  nv 
(CFFS), 
fish  nr3 

No.  of 

independent 
samples  (n) 

Least  squares 

TSqawrs 

(dB  re  1  m2  at  415  Hz) 

14 

12:17 

3000 

0.04 

85 

-39.7 

14 

14:03 

3500 

0.03 

59 

-40.1 

15 

10:08 

210 

0.02 

12 

-40.3 

15 

11:19 

105 

0.02 

6 

-40.3 

15 

12:20 

87 

0.02 

2 

-39.8 

15 

13:19 

306 

0.02 

17 

-39.8 

Mean  TSqawrs 

-40.0 

SD  TSqawrs 

0.7 
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Appendix  F.  Signal  to  noise  ratio  in  OAWRS  intensity  images 


The  signal  to  noise  ratio  (SNR)  in  OAWRS  images  is  the 
ratio  of  the  expected  scattered  intensity  from  fish  to  that 
from  the  seafloor 


SNR  =  101og10 


Ofos^fish) 

bottom) 


(FI) 


Expressions  for  the  second  moment  of  the  scattered  field 
from  fish  groups ,  ( 1 4>s  1 2  )fish ,  proportional  to  the  expected  scat¬ 
tered  intensity,  appear  in  Appendix  C. 

The  second  moment  of  the  scattered  field  from  the 
seabottom,  (l^sl2)bottom  r  depends  on  seafloor  scattering  prop¬ 
erties.  Seafloor  scattering  has  been  studied  extensively 


since  World  War  II,  and  empirical  models  have  been  used  to 
describe  various  field  measurements  (Urick  1983).  For  most 
seabottom  types,  including  sand,  silt  and  rocky  bottoms,  lit¬ 
tle  or  no  frequency  dependence  has  been  measured  in 
seafloor  scattering  strength  over  the  OAWRS  frequencies 
considered  here  (Urick  1983).  Seafloor  scattering  properties 
measured  recently  (Galinde  et  al.  2008)  are  used  for  all  con¬ 
tinental  shelves  considered  here  for  a  typical  sandy  bottom. 
The  seafloor  returns  for  all  the  environments  discussed  in 
'Potential  ecosystem  exploration'  in  the  main  paper  are 
computed  by  using  the  scattering  model  developed  in 
Galinde  et  al.  (2008). 
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